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We present new results from our test of Lorentz invariance, which compares two orthogonal 
cryogenic sapphire microwave oscillators rotating in the lab. We have now acquired over 1 year 
of data, allowing us to avoid the short data set approximation (less than 1 year) that assumes 
no cancelation occurs between the ke- and ko+ parameters from the photon sector of the standard 
model extension. Thus, we are able to place independent limits on all eight Ke_ and ko+ parameters. 
Our results represents up to a factor of 10 improvement over previous non rotating measurements 
(which independently constrained 7 parameters), and is a slight improvement (except for ftfj?) over 
results from previous rotating experiments that assumed the short data set approximation. Also, an 
analysis in the Robertson-Mansouri-Sexl framework allows us to place a new limit on the isotropy 
parameter Pmm = i5 — /9 + ^ of 9.4(8.1) x 10~^^, an improvement of a factor of 2. 

PACS numbers: Valid PACS appear here 



In recent times there has been an increase in activity 
in experimental tests of Local Lorentz Invariance (LLI), 
in particular light speed isotropy tests with at least 6 
experiments reported in the last 3 years P, S S IE ■ 
This is largely due to advances in technology, allowing 
more precise measurements, and the emergence of the 
Standard Model Extension (SME) as a framework for the 
analysis of experiments, providing new interpretations of 
LLI tests. None of these experiments have yet reported 
a violation of LLI, though the constraints on a putative 
violation have become more stringent by approximately 
three orders of magnitude in the same time frame. 

LLI is an underlying principle of relativity, postulating 
that the outcome of a local experiment is independent of 
the velocity and orientation of the apparatus. Tests of 
LLI are motivated by the central importance of this pos- 
tulate to modern physics, as well as the development of a 
number of conflicting unification theories, which suggest 
a violation of LLI at some level. To identify a viola- 
tion it is necessary to have an alternative theory to inter- 
pret the ex p eriment [7 1 , and many have been developed 
llliliSllljllElCTIl^. The kinematical frameworks 
(RMS) [Sl ISl postulate a simple parametrization of the 
Lorentz transformations with experiments setting limits 
on the deviation of those parameters from their values 
in special relativity (SR). Due to their simplicity they 
have been widely used to interpret many experiments 
[E IM IE 0| . More recently, a general Lorentz vio- 
lating extension of the standard model of particle physics 
(SME) has been developed [Ullllli whose Lagrangian 
includes all parameterized Lorentz violating terms that 
can be formed from known fields. This has inspired a 
new wave of experiments designed to explore uncharted 
regions of the SME Lorentz violating parameter space. 

Our experiment consists of two cylindrical sapphire 



resonators of 3 cm diameter and height supporte d by 
spindles within super-conducting niobium cavities |18l |. 
The sapphire loaded cavities are situated one above the 
other, oriented with their cylindrical axes orthogonal to 
each other in the horizontal plane. The experiment is 
rotated with a period of 18 seconds around its vertical 
axis. Whispering gallery modes [li| are excited in each 
near 10 GHz, with a difference frequency between the 
two of 226 kHz. The difference frequency along with var- 
ious experimental parameters are logged by a stationary 
data acquisition system as a function of the experiments 
orientation. A detailed description of the experiment can 
be found in 20]. 

Inside the sapphire crystals standing waves are set up 
with the dominant electric and magnetic fields in the ax- 
ial and radial directions respectively, corresponding to 
a Poynting vector around the circumference. The fre- 
quency of each resonator v is proportional to the speed 
of light c and inversely proportional to the electrical path 
length L of the resonator [v oc c/L), where L is depen- 
dent on the material properties of the sapphire crystal, 
which have been shown to have a negligible dependence 
on orientation 21] . Hence, by measuring the difference 
frequency between the two orthogonal cavities as they 
rotate we make a direct observation of the isotropy of 
the speed of light. 

To test for Lorentz violations we derive the perturba- 
tion of the difference frequency with respect to an al- 
ternative test theory. In the photon sector of the SME 
this may be calculated to first order as the integral over 
the non-perturbed fields (Eq.(34) of HE , see [ll for 
an application to our case). The change in orientation 
of the fields due to the rotation of the experiment in 
the lab and Earth's orbital and sidereal motion induces 
a time varying modulation of the difference frequency. 
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which is searched for in the experiment. In the photon 
sector of the SME Lorentz violating terms are pa- 
rameterized by 19 independent components, which are in 
general grouped into three traceless and symmetric 3x3 
matrices (Ke+, ^o-, and Ke-), one antisymmetric matrix 
(ko+) and one additional scalar, which all vanish when 
LLI is satisfied. The 10 independent components of Ke+ 
and Ko- have been constrained by astronomical measure- 
ments to < 2 X 10^'^^ 0,1131 • Recently two combina- 
tions of these parameters have been further constrained 
to less than parts in lO"*^^ |2^ . The scalar ktr component 
has been constrained to < 10"'' by 24| through the re- 
analysis of previous Ives-Stilwell experiments, who also 
propose interferometric techniques to improve on this by 
seven orders of magnitude. Seven components of 
and Ko+ have been independently constrained in station- 
ary optical and microwave cavity experiments |1|, iSj at 
the lO^'^^ and lO^'^^ level respectively. The last remain- 
ing component k^j^ was only recently constrained for the 
first time by a group of cavity experiments [SJi[Jp.j25L 
designed to both improve on the results of [l|, |2, |3| and, 
more importantly, be sensitive to k^_^ through the use of 
active rotation in the laboratory. 

However, the most stringent independent limits on 
the isotropy (ke-) and boost terms (ko+) can only be 
achieved with 1 year of data. This is because the maxi- 
mum boost with respect to the Sun Centered Equatorial 
Celestial Frame (SCECF) is due to the Earth's annual 
motion. Thus, over 1 year of data is required to decor- 
relate the parameters. Previous analysis which 
contained significantly less than one year of data, con- 
strained the Kg- and Ko+ parameters by assuming no 
cancelation occurred in the case of a non-zero Lorentz 
violating effect. We have now acquired sufficient data to 
remove this assumption, producing independent limits on 
all of the eight components of Kg- and Ko+- 

Alternatively, with respect to the RMS framework, we 
analyze the change in resonator frequency as a function 
of the Poynting vector direction with respect to the ve- 
locity of the lab in some preferred frame (as in [l7ll20| ). 
typically chosen to be the cosmic microwave background. 
The RMS parameterizes a possible Lorentz violation by 
a deviation of the parameters (a, /3, S) from their SR val- 
ues (—^,^,0). These are typically grouped into three 
linear combinations representing a measurement of (i) 
the isotropy of the speed of light {Pm m = S — P + i), 
a Michelson-Morley (MM) experiment [23, constrained 
by to parts in 10^^ (u) the boost dependence of the 
speed of light {Pkt — (3 — a — 1), a, Kennedy-Thorndike 
(KT) experiment ||22|], constrained by 0| to parts in 
10"'', and (iii) the time dilation parameter (P/5 = ct +^) , 
an Ives-Stillwell (IS) experiment [2^ , constrained by |30 | 
to parts in 10"^. Because our experiment compares two 
cavities it is only sensitive to Pmm- 

In our previous analysis the amplitude and phase of 
a Lorentz violating signal was determined by fitting the 



parameters of Eq^to the data, with the phase of the fit 
adjusted according to the test theory used. 

= A+Bt+S^ d cos{uJit+ipi)+Si s\n{LjJit+(pi) (1) 

Here is the average unperturbed frequency of the 
two sapphire resonators, and Ai/q is the perturbation of 
the 226 kHz difference frequency. A and B determine the 
frequency offset and drift, and Ci and Si are the ampli- 
tudes of a cosine and sine at frequency uji respectively. 
In the final analysis we fit 15 frequencies to the data, 
LUi = {2ujR, 2uiR ± VL,^,2ujr ± uit^,2u!R ± ± VLi^,2uir ± 
2[jj0, 2uj R±2uj(^±^l@) , where wr is the rotation frequency 
of the experiment in the lab and and Jl® are the side- 
real and annual frequencies of the Earth's rotational and 
orbital motion respectively. Since the residuals of the 
fit exhibit a significantly non- white behavior, the opti- 
mal regression method is weighted least squares (WLS) 
0|- WLS involves pre- multiplying both the experimen- 
tal data and the model matrix by a whitening matrix 
determined by the noise type of the residuals of an or- 
dinary least squares analysis. However, this method of 
analysis proved to be computationally intensive due to 
the large amount of data we have now acquired. For 
this reason, an alternative approach used by 0, @ was 
adopted. Using this technique we reduce the size of the 
data set by demodulating it in quadrature with respect 
to 2u)R in blocks of 40 periods of rotation. The num- 
ber of periods was chosen to minimize the net effect of 
narrow band noise (due to instabilities in the system- 
atic at 2u!r) and broad band noise (due to oscillator fre- 
quency noise), which is similar to an optimal filter. By 
fitting the expression of Eq[21to each block of data using 
an ordinary least squares regression technique we deter- 
mine the coefficients S{t) and C(i), which can be con- 
sidered linear combinations of the sidereal, semi-sidereal, 
and annual modulations and combinations thereof. The 
relationship between S{t) and C{t) and the various mod- 
ulation frequencies is expressed in EqsOl and ^ where 
LOi = (ri0, ± VIq,2lu^, 2u)^ ± il®). 



A + Bt + S{t) siTL{2ujRt + + C{t) cos{2LURt + (p) 

(2) 



S{t) = 5o +^ S's.i sm{ujit + (pi) + Sc,i cos{ujit + ipi) (3) 



C{t) = Co+'^Cs,i sm{uJit + (pi) + Cc,z cos{uJit + ipi) (4) 

i 

A comparison was made between the two techniques 
by performing a complete analysis of 30 data sets (3 data 
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FIG. 1: Cosine, C{t), and sine, S{t), amplitudes resulting 
from demodulation of the data at 2ujr in blocks of 40 rota- 
tions, with a linear fit removed from each data set. 
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sets were later excluded from the analysis due to overly 
large and varying systematic signals at 2wj^). Both tech- 
niques produced consistent results, with the uncertainties 
associated with the demodulated technique being lower 
than the WLS technique by no more than 15 percent. 
The difference between the two techniques is most likely 
due to the efficiency with which the data analysis could 
be optimized for the noise type present in the data. WLS 
only takes into account the broad band noise (spectral 
density) whereas the optimization used in the demodu- 
lated technique takes into account the extra noise source 
of instability of the systematic at 2ujfi. Hence, the lat- 
ter approach was adopted in further investigations of the 
data. 

The data used in this analysis spans a period from De- 
cember 2004 to January 2006. It consists of 27 sets of 
data totalling approximately 121 days. Shown in Fig. 
n are the S{t) and C{t) resulting from the demodula- 
tion of the data at 20;^. An offset and drift has been 
removed from the coefficients derived from each data 
set. As described earlier, this data is then used to de- 
termine the amplitudes of the frequencies of interest. In 
we describe how systematic effects dominate the data 
at 2w/{, limiting our ability to constrain test theory pa- 
rameters associated with this frequency (a detailed dis- 
cussion of the systematics and their effect is thus left out 
here). Also, we do not consider the nearby annual offsets 
(2a;ij ± fi©) for two reasons. Firstly, the strong system- 
atic signal at 2ujr has been shown to have a significant 
effect on nearby sidebands due to leakage 0, and sec- 
ondly, by subtracting a linear drift from the individual 
data sets after being demodulated (as presented in FigCJ 
it is possible that a signal at the annual frequency may 
be suppressed. However, all other frequencies of interest 
(see TablD are close to the sidereal or semi-sidereal fre- 
quencies, so will be unaffected by the removal of an offset 



FIG. 2; The amplitudes So and Co for each of the 27 data 
sets (squares), used to limit the parameter ftfj? of the SME. 
Also shown (circle) is the mean and its standard error (5*0 = 
2.9(4.6), Co = -3.0(3.8) x 10"^^). 



and drift from each data set. 

In the SME, all ke- and K0+ parameters other than 
k^- can be constrained from the sidereal and semi- 
sidereal frequencies and their annual frequency offsets as 
outlined in Tab|l| k^^ only appears in the coefficient 
C2ujrt SO to determine a limit we do the same as in |^ 
and consider the C2ujii coefficients for each data set to 
be independent and treat them statistically. The sys- 
tematic at 2ujii has been shown to be primarily due to 
tilt variations. It remains relatively constant in phase 
within a data set but varies between data sets. Fig. |21 
shows the C2ujr and S2ujr coefficients for the 27 data sets. 
Also shown is the mean and standard error of the mean 
which is used to calculate k^^ . The results for the SME 
analysis are given in Tab ITU We note that the results 
for k^^ and k-^_^ are significant at approximately the 
3(7 and 2a level respectively. However, we do not be- 
lieve this to be an indication of a Lorentz violating effect 
for reasons similar to those given in 0], which also used 
data taken over more than one year. Our result for k^_^ 
is inconsistent with other recent measurements shown in 
Tabim Also, an examination of the corresponding side- 
band coefficients from an analysis of the individual data 
sets (not shown here) shows no coherence in the phase of 
the signal, which would be expected in the presence of a 
genuine Lorentz violating effect. 

In terms of the RMS framework, the advantage to be 
gained by having one year of data is primarily statis- 
tical. Due to the symmetry of our experiment we are 
not sensitive to the boost parameter of the RMS, Pkt, 
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TABLE I; Shown are the relationships between the ke- and Ko+ parameters of the SME and the coefficients Cc,uji,Cs,uii,Sc,u}i 
and Ss,uii from Eqs. |31and|l]for the 8 frequencies of interest, normalized for the experimental sensitivity S. x is the colattitude 
of the lab, rj is the declination of the Earth's orbit relative to the SCECF, and /3® and /3l are the boost suppression terms 
due to the lab velocity from Earth's orbital and rotational motion respectively. Also shown is the measured value (in 10~^^) of 
each coefficient used in the analysis along with its statistical uncertainty. The values for Cc,o (used to constrain wfj?) and Sc,o 
were determined by averaging over the data sets (see text). The coefficients of fJ® were not included in the analysis (see text). 



fsin^(x)Sfi -30(38) 
-/3®sin2(x)(cos(r?)Kf/ + 2sm{r])k^^) 

/3®cos(x)sin(x)sin(r?)Kr-f -2.3(0.7) 

sin(2x)ftf_^ + 2/3z,Kf/ 1.9(0.7) 

/3®cos(x)sin(x)sin(»?)Srf -2.0(0.7) 

-i/3®cos''f (3 + coai2x))ko+ -0.4(0.7) 

-i(3 + cos(2x))(Sfi^ - kD -0.6(0.7) 

^/3ffisin'a(3-Hcos(2x))af/ -3.4(0.7) 



-/3esin2(x)Kf/ 
-2/3g5Cos(f )sin(x)(cos(f )Sf/ - sin(f )Sf/) 

sm{2x)kl^ +2f3Lk^f 
-/3esin(f)sin(2x)(cos(f)K^^ 



-"i/3®cos2f(3-^cos(2x))srf 
■i(3 + cos(2x))Kf^ 



+ smm)kf/) 



ffism^ f (3 - 



' cos(2x))fto 



YZ 



0.9(0.7) 
-2.5(0.7) 
-1(0.7) 
-0.7(0.7) 
-1.7(0.7) 
-0.5(0.7) 



/3ecos(2)sin(2x)(cos(f )Aif/ - sin(2)sf/) 

-2(sin(x)Kr- + 2/3iCOs(x)Kr+^) 
2/3esin(2)sin(x)(cos(2);if/+sin(f)^f/) 

2/3ecos2|cos(x)Sr-f 

2cos(x)KfJ' 
-2/3esin2fcos(x)Srf 



29(46) 

-0.3(0.8) 
1.4(0.8) 
-5.4(0.8) 
1.5(0.8) 
-1.2(0.8) 
0.5(0.8) 



/3®sin(x)sin(r;)K^/ 
2(sin(x)Kf/ + 2f3LCOs{x)k^+) 
/3esin(x)sin(r;)K^/ 
-2/3ecos2|cos(x)Kf/ 
-cos{x){k^J' -klJ) 
2/3ffiSin2 2cos(x)Aif/ 



-0.8(0.8) 
-3.6(0.8) 
0.4(0.8) 
-3.2(0.8) 
2.8(0.8) 
3.4(0.8) 



and cavity experiments are not sensitive to the time di- 
lation parameter a. However, we can improve on our 
previous constraint on the isotropy parameter Pmm by 
taking a weighted average over the results of multiple 
data sets. We analyze each data set using the WLS tech- 
nique described earlier. The association between Pmm 
and the coefficients of the frequencies of interest is de- 
scribed in 5j. The coefficients of Eq^are for the frequen- 
cies oJi = {2ljJr, 2ujr±ujq, 2u!ii±2uji^) only. We calculate 
a value for the RMS parameter of 9.4(8.1) x 10~". 

In conclusion, by collecting over one year of data we 
have been able to set the first independent limits on 8 
parameters in the photon sector of the SME, without as- 
suming that no cancelation occurs between the isotropy 
and boost terms. The results do not indicate any Lorentz 
violating effects, and compared to previous experiments 
we see a slight improvement in the constraints on these 
parameters. We improve on our previous determination 
of Kf_5 by more than a factor of three. However, due to 
the systematic disturbances present at twice the rotation 
frequency we are unable to measure this parameter with 
the precision of 6] , who have developed a tilt control sys- 
tem which avoids the major rotation induced systematic. 
Also, we have reduced the limit on the isotropy parame- 
ter Pmm of the RMS framework by a factor of two. 

To improve on these results we intend to replace the 
resonators with higher quality sapphire loaded cavities, 
which have a frequency instability approximately 40 
times lower than the current experiment [3l| . Consid- 
erable effort will need to be invested to improve the ro- 



TABLE 11: Results for the SME Lorentz violation parame- 
ters determined independently in this work. Also shown for 
comparison is the previous best independent constraints of 
seven parameters Q and more recent short term results that 
assume no cancelation between the Ke- and ko+ terms, other 
than Kff |l@(Ke- in 10"^^ ko+ in 10"^^). The Pmm pa- 
rameter from the RMS framework is also listed (in 10~^^). 
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2.1(2.1) 


-5(13) 


-4.5(3.7) 


{kg- ~ '^e- ) 
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k^F 


-4.4(2.5) 


-14(23) 


-3.6(2.7) 


k^^ 


-3.2(2.3) 


27(22) 


2.9(2.8) 


Pmm 


9.4(8.1) 


120(220) fl7] 


-21(19) 



tation system and reduce environmental disturbances for 
this improvement to be realized. 
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